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Abstract 

We theoretically investigate a non-magnetic impurity effect on the temperature dependence of the vortex core shrinkage 
(Kramer-Pesch effect) in a single-band s-wave superconductor. The Born limit and the unitary limit scattering are 
compared within the framework of the quasiclassical theory of superconductivity. We find that the impurity effect 
inside a vortex core in the unitary limit is weaker than in the Born one when a system is in the moderately clean 
regime, which results in a stronger core shrinkage in the unitary limit than in the Born one. 
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1. Introduction 

The radius of a vortex core in type-II supercon- 
ductors is one of the fundamental physical quantities 
which characterize a property of superconductivity. 
The temperature and magnetic field dependence of 
the core radius has been investigated theoretically 
and experimentally [1-19]. The low-temperature 
vortex core shrinkage, called the Kramer-Pesch 
(KP) effect [1], was theoretically investigated un- 
der the influence of non-magnetic impurities in the 
Born limit previously [2] . Impurity effects are char- 
acterized by the scattering phase shift related to the 
impurity potential strength [20-23]. The Born limit 
corresponds to the limit of weak impurity potential 
and correspondingly small phase shift. The opposite 
limit is called the unitary limit, where the impurity 
potential is infinitely strong and the phase shift 
is 7r/2. The difference between these limits plays 
an important role in, for example, unconventional 
superconductors [20-22] . 

In this paper, we theoretically study the KP ef- 



fect both in the Born and the unitary limit in an 
s-wave superconductor, and compare their results. 
It is found that the temperature dependence of the 
core shrinkage is stronger in the unitary limit than 
in the Born one, in the moderately clean regime 
wliere the mean free path is of the order of or larger 
than the coherence length. Such a difference of the 
core shrinkage can be investigated experimentally 
by, e.g., muon spin rotation [3,5,6], scanning tunnel- 
ing microscopy [24] , resistivity [4] , and specific heat 
[11,25] if there is a suitable superconducting mate- 
rial in which different types and densities of impu- 
rities can be doped. 

2. Formulation 

We consider a single vortex in a single-band s- 
wave superconductor. The system is assumed to be 
an isotropic two-dimensional conduction layer per- 
pendicular to the vorticity along the z axis. In a cir- 
cular coordinate system within the layer, the real- 
space position is r = (r cos 0, r sin 0). The unit vec- 
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Fig. 1. Schematic figure of the pair potential as a function 
of the distance from the vortex center. It depicts the length 
§1 that characterizes the vortex core radius [see Eq. (1)]. 

tor fe represents the sense of the wave number on 
a Fermi surface assumed to be circular. The Fermi 
velocity is Vp — v^k. The pair potential around 
the vortex is A(r) = A{r,(p) ~ |A(r)| exp(i(/)). We 
will consider the temperature T dependence of the 
length ^1 that characterizes the vortex core radius 
[1-3], 
1 
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This quantity is depicted in Fig. 1. Note that is 
related to the pair potential slope at the vortex cen- 
ter and scales with the distance at which the vortex 
current reaches its maximum value [1,3,6,11], while 
|A(r)| is restored at a distance of the order of the 
coherence length {'S> £,i in the clean limit) even at 
low temperatures [11,26]. 

To obtain (T) , the vortex core structure is cal- 
culated by means of the quasiclassical theory of 
superconductivity as in Ref. [2]. The Eilenberger 
equation is numerically solved to obtain the qua- 
siclassical Green's function g(iu)n,r,k). The effect 
of impurities distributed randomly in the system is 
taken into account through the impurity self energy 
S(ia;„,r,fc). The quasiclassical Green's function, 
the impurity self energy, and the Eilenberger equa- 
tion to be solved are, respectively, given as [2,27-31] 
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(3) 



ivp ■ Vg + [ioJnTs ~ \g] = 0. 

The equation is supplemented by the normalization 
condition g^ ~ — tt^to [29,32]. Here, fa is the z com- 
ponent of the Pauli matrix, tq is the unit matrix. 
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Fig. 2. Spatial profiles of the pair potential amplitude |A(r)| 
around the vortex under the influence of impurity scattering 
in the unitary limit. The horizontal axis r is the distance from 
the vortex center. The scattering rate is Fn/Ao = 0.1 (solid 
lines) and Fn/Ao = 1 (dashed lines). For each scattering 
rate, the temperature is T/Tc = 0.1-0.7 from top to bottom 
by 0.2 step. 

and the brackets denote the commutator [A, B] = 
AB — BA. The Eilenberger equation contains the 
renormalized Matsubara frequency (pair potential) 

cj„ (A) defined by 



A = 





-A* 



Sd, 

A 





(4) 
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We consider an isolated single vortex in an extreme 
type-II superconductor (Ginzburg-Landau parame- 
ter K ^ 1), and therefore the vector potential is 
neglected in Eq. (3). Throughout the paper, we use 
units in which h = kB — I- 

The Eilenberger equation (3) can be solved by the 
Riccati parametrization [33-35]. The quasiclassical 
Green's function is expressed as 
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I — ab i2a 
-i2b -{1-ab) 



(6) 



Here, sgn(aj„) is the signum (or sign) function. The 
two quantities aiiuJn, r, k) and 6(iw„, r, fe) are inde- 
pendently determined by solving the Riccati equa- 
tions, 



vf ■Va+ {2ojn + A*a)a - A = 0, 
vp-Vb- {2(:jn + Ab)b + A* = 0. 



(7) 
(8) 
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Fig. 3. (a) The temperature T dependence of the vortex core 
radius ^i. (b) The scattering rate r„ dependence of at 
T/Tc = 0.02. Note that the horizontal axis is in log scale in 
(b). The lines are guides for eyes. 

These differential equations are solved along a 
straight line parallel to vp [33,36,37] by using the 
bulk solutions as initial values [33,38], 



b+oo — ■ 
and 

a+oo = 



Z + |A| 



(w„ > 0), 



(Wn > 0), 



(9) 
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Z + |A| 
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A* 



c^2 + |A|- 



K <0),(11) 



K < 0).(12) 



A stable numerical solution for a {b) is obtained by 
solving the Riccati equation in forward (backward) 
direction along the straight line for a;„ > [33,39]. 



By contrast, the equation for a (b) is solved in back- 
ward (forward) direction for w„ < 0. 

Considering an s-wave non-magnetic impurity 
scattering and the i-matrix, S is given by [2,27,40] 



T,{iuj„,r) 



r„ 



l-(sin^<5o)(l-C) 
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where C = (5)^ -I- {f){P) with (• • •) being the av- 
erage over the Fermi surface with respect to k. The 
impurity scattering rate in the normal state is Fn, 
which is related to the mean free path I — UF/2Fn. 
The scattering phase shift is Sq- We set (5o = in the 
Born limit (keeping Fn finite) and Sq = 7r/2 in the 
unitary limit. 

The self-consistency equation for A, called the gap 
equation, is given as 



A(r) = XttT 



E 



(14) 



where Wc is the cutoff energy and the coupling con- 
stant A is given by 



1 
A 



= In 



E 



0<n<(wc/7rT-l)/2 
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(15) 



Here, Tc is the superconducting critical tempera- 
ture. We set uic = lOAo with Aq being the BCS 
pair-potential amplitude at zero temperature. 

The Eilenberger (Riccati) equation, the impurity 
self energy, and the gap equation are numerically 
solved self-consistently. The used boundary condi- 
tions for the pair potential and impurity self energy 
far from the vortex are the same as those discussed 
in Ref. [2]. See the Appendix for more details on 
the calculation procedure. In the next section, we 
will show results obtained from self-consistent so- 
lutions. We define the zero-temperature coherence 
length ^0 = vf/Aq. 

3. Result and Discussion 

We show numerical results for the pair potential 
A(r) in Fig. 2, which are obtained under the in- 
fiuence of impurity scattering in the unitary limit. 
Similar results are obtained in the Born limit (not 
shown). The vortex core radius ^1 is related to the 
slope of A(r') at the vortex center r — according to 
Eq. (1) (see also Fig. 1). As seen in Fig. 2, the slope 
increases with decreasing the temperature T. The 
slope becomes smaller with increasing the scattering 
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(a) Born limit (T/Tj,=0.02, rr,/Ao=0.1) (a) Born limit (1/1^=0.02, rn/Ao=0.1) 




Fig. 4. Spatial profiles of the imaginary part of tfie diagonal 
impurity self energy for several Matsubara frequencies 
LUn in the Born limit (a) and in the unitary limit (b). The 
horizontal axis r is the distance from the vortex center. Note 
that ^i/5o ~ 0.13 in the Born limit and 0.05 in the unitary 
one at T/Tc = 0.02 and r„/Ao = 0.1. 

rate r„ at the low temperature T/Tc = 0.1, which 
indicates the suppression of the vortex core shrink- 
age due to the impurity scattering. We investigate 
the impurity effect on the vortex core radius below. 

The temperature dependence of the vortex core 
radius [Eq. (1)] is shown in Fig. 3(a). £,i{T) is 
plotted for the scattering rate r„/Ao = 0.1 and 1 in 
the Born and unitary limits. The KP effect means 
^1 cx T and ^ in T — ^ [1]. The suppression 
of the KP effect, namely the weakened vortex core 
shrinkage due to impurities, is seen in Fig. 3(a) . That 
is, with increasing the scattering rate r„, the slope 
of £,i{T) becomes small and the intercept value of 
^1 in the limit T becomes large. Exceptionally, 
the plot for r„/Ao = 0.1 in the unitary limit seems 
to exhibit a rather strong KP effect. To see it in 
more detail, in Fig. 3(b) we plot the dependence of 
^1 on r„ in the Born and unitary limits at T/Tc = 



Fig. 5. Spatial profiles of the real part of the off-diagonal 
impurity self energy S12 for several u!„ in the Born limit 
(a) and in the unitary limit (b). r is the distance from the 
vortex center along the radial line in the direction </< = 0. 

0.02, which is the lowest temperature in Fig. 3(a). 
The difference between the Born and unitary limits 
is clearly seen in the moderately clean regime (l > 
^0 > ^Ao/ep, namely Ao/2£f < r„/Ao < 0.5;£Fis 
the Fermi energy). This result suggests that inside a 
vortex core the impurity effect is less effective in the 
unitary limit than in the Born one when the system 
is neither dirty nor too much clean. 

As a clue to the difference between the Born and 
unitary limits, in Figs. 4 and 5 we show the impu- 
rity self energies and S12 as functions of the dis- 
tance r from the vortex center. Here, we plot ImE^ 
and ReEi2 [41] for, as representatives, several en- 
ergies uJn — TrT(2n + 1) (n — 0, 1, 2) and ^139 = 
2797rT ~ Wc at T/T^ = 0.02. The scattering rate is 
set r„/Ao = 0.1 in common. We notice that the im- 
purity self energies are enhanced in the vicinity of 
the vortex center for small |w„| in the Born limit, 
but not in the unitary one. 

Then, let us investigate the relation between the 
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behavior of the self energies and the impurity effect. 
Impurity effects are ineffective under the condition 
that the impurity self energies Sd and — S12 are pro- 
portional to iujn and A, respectively, with a com- 
mon proportionality constant [20,42]. It is the so- 
called Anderson's theorem [43]. Here, we introduce 
a dimensionless quantity that measures the degree 
of the deviation from the above condition: 



(a) Born limit (T/T(;=0.02, rn/Ao=0.1) 



6a ^ 



A(r) 



(16) 



We show Sa in Fig. 6, where r„ are set the same 
between the Born and unitary limits. Indeed, Sa is 
zero far away from the vortex center, where Ander- 
son's theorem is satisfied. In contrast, Sa becomes 
finite inside the vortex core. It is clearly seen that 
i5a is prominently large near the vortex center in the 
Born limit, compared with the unitary-limit case. 
This result indicates that the impurity effect inside 
the vortex core is weaker in the unitary limit than 
in the Born one, resulting in the weaker suppression 
of the KP effect in the unitary limit (Fig. 3). 

4. Conclusion 

We studied the non-magnetic impurity effect on 
the low-T vortex core shrinkage in the s-wave su- 
perconductor. We found that in the moderately 
clean regime, the suppression of the core shrinkage 
is weaker in the unitary limit than in the Born one. 
It is attributed to the difference in the impurity 
self energy between the Born and unitary limits. 
However, it is difficult to intuitively understand the 
reason why the difference appears. Instead, we in- 
troduced the indicator Sa to estimate the efHciency 
of the impurity effect. From the analysis of Sa, it 
was elucidated that the impurity effect inside the s- 
wave vortex core is stronger in the Born limit than 
in the unitary one in the moderately clean regime 
at low temperatures [44] . 

Appendix 

Here, we briefly explain the calculation procedure 
on the single vortex. Our starting point is the Eilen- 
berger equation (3) with the impurity self energy 
and the impurity-averaged Green's function (refer 
to Refs. [29-32] for its derivation). It has been used 
to investigate spatially inhomogeneous systems such 
as vortices [2,17-19,27,28,31,45]. The Eilenberger 



< 

to 




(b) unitary limit (T/Tj,=0.02, rn/Ao=0.1) 




Fig. 6. Plots of the indicator 5a [Eq. (16)] as a function of 
the distance r from the vortex center for several ujn in the 
Born limit (a) and in the unitary limit (b). The scattering 
rate is set r„/Ao =0.1 in common so as to compare the 
data between (a) and (b). 

equation is transformed to the Riccati equations (7) 
and (8) (see Refs. [33,35] for the details). The Eilen- 
berger equation corresponds to three coupled differ- 
ential equations [46] , and therefore it is not so easy 
to solve. It is easier to solve the Riccati equations, 
which are two decoupled differential equations. Sta- 
ble solutions are easily obtained when solving the 
Riccati equations numerically [39]. 

The Eilenberger and the Riccati equation contain 
the differential operator in the form Vp ■ 'V, which 
is treated as follows. Consider the orthogonal coor- 
dinate system r = xx + yy +zz with z parallel 
to a rectilinear vortex line, {x, y, z are orthogonal 
unit vectors.) Here, the vortex center is situated on 
the z axis. We express the Fermi velocity as vf = 
vpk — VF{cos(j)k s'mOkX + sin(f)k siiiBky + coaOkz). 
Because of the translational symmetry along the 
vortex line, one can omit the z dependence. Thus, 
■ V ^ Vyi_ ■ V = xvYC0S(f)kSin9k{d/dx) + 
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y vf sin (p k sin 9 k{d/dy), where vf± is the Fermi ve- 
locity projected onto the xy plane normal to the vor- 
tex line. When the Fermi surface is a cylinder with 
central axis parallel to the z axis (|| the vortex line), 
the angle 9k is set to tt/2. It is convenient to in- 
troduce a new orthogonal coordinate system where 
one of the axes is parallel to Vf±. Using such coor- 
dinates {s,t), r = XX + yy = ss + tt with s \\ vf±. 
The variable t is called the impact parameter. Even- 
tually, t)F • V = sin 9k{d/ds) in the equations to 
be solved. One can obtain the quasiclassical Green's 
function at a position rg = xqX + yoy = sqS + 
tot by solving the first-order differential equations 
along a straight line parallel to the s axis with fixing 
t = tf). The k dependence is obtained by solving the 
equations in the same way for each different direc- 
tion s and polar angle 9k ■ Note that s is a function 
of the azimuth angle <pk of fc. 

A single vortex is characterized by the pair po- 
tential A(r) = |A(r)| exp(i0) in the case of an 
s-wave superconductor. Because of the rotational 
symmetry around the vortex line, the dependence 
on the real-space azimuth angle </> is factored out 
as / = /exp(i(/)), p =_ p exp{-i(j)), g = g, 
Si2 = Si2exp(z0), i;2i_= S21 exp(-i(/)), Ed = Ed, 
a = aexp(i0), and b = bexp{—i(j)), where the quan- 
tities with "bar" are independent of (j). Therefore, 
the task is to determine the radial r dependence. 

In numerical calculations, the pair potential and 
the impurity self energies are represented by their 
values at discrete points on the radial line. The val- 
ues between those positions are calculated by linear 
interpolation. It is necessary to take dense point 
spacing near the vortex ccinter because the values 
have rapid spatial variations there, while sparse 
points are sufficient far from a vortex. In consid- 
eration of it, we take non-equally-spaced discrete 
points represented by r, = f[exp(i?j) — l] with 
equally-spaced discrete points Ri in dimension-less 
space. In the present study, we set f = 10~^^o and 
take 41 points from the vortex center to the cutoff 
distance Vc = lO^o- In the outside region r > Vc, 
we set |A(r > rc)| |A(r = rc)|, Si2(r > r^) = 
Si2(r = rc),_S2i(r > Tc) = T,2i{r = Tc), and 
Sd(r > re) = Sd(r = re) [2]. 

The equations are numerically solved itcra- 
tively until the self consistency is attained as 
< 2 X 10~^, where 6Q = max^j |(3new(ri) — 
Qoid(ri)|/|<5new(ri)|} and Q stands for the pair po- 
tential and the self energies. The larger the scatter- 
ing rate Fn is, the slower the convergence becomes. 



especially in the unitary limit. Therefore, we use 
an acceleration method [47] for updating the pair 
potential and the self energies at each iterative step. 
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